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ABSTRACT 

We use a recent sample of 49 galaxies to show that there is a proportionality relation between 
the black hole mass M Bii and the quantity ji = Mqcj/c, where Mq is mass of the spheroidal 
stellar component and a is the stellar velocity dispersion, ji is called the momentum parameter 
and the ratio is Mbh/jU ~ 3.3. This result is applied to the penetrating-jet feedback model 
which argues that the correlation that holds is with a momentum-like parameter, although this 
feedback mechanism is based on energy balance. 
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1 INTRODUCTION 



Relations between a supermassive black hole (SMBH) mass, 
Mbh, and other properties of its host galaxy have been stud- 
ied intensively in the last decade. Two of the galactic properties 
most commonly correlated with Mbh are the stellar mass of the 
spheroidal component (which we refer to as the bulge), Mq (e.g. 
Kormendy & Richstone 1995; Magorrian et al. 1998; Laor 2001; 
Hu 2009; Graham & Spitler 2009), and it's stellar velocity dis- 
persion, o~ (e.g. Gebhardt et al. 2000; Merritt & Ferrarese 2001; 
Graham 2008a; Graham (2008b); Hu 2008; Shen et al. 2008; 
Giiltekin et al. 2009). These relations are often assumed to have the 
form of a power law (i.e. linear when plotted on a log-log scale). 
However, due to the large scatter there is still no consensus on 
the best parameters of the different relations. For example, despite 
some claims for a proportionality relation between the SMBH and 
bulge mass, Laor (2001) found that the ratio Mbh/Mq increases 
with mass. This conclusion was strengthened by some models and 
simulations (e.g. Shabala & Alexander 2009). 

The more significant disagreement is on which fundamen- 
tal galactic property is behind the correlations. In a recent paper, 
Feoli et al. (2010) studied three samples of SMBH masses and their 
host galaxies. They argued that Mbh is better correlated with the 
energy parameter, MqO 2 , than with Mq or with o~ alone. Behind 
this comparison stands the view that the feedback between the ac- 
creting SMBH and its environment is driven by energy. 

Results from recent years show that the process of galaxy for- 
mation requires another energy source to that of the gravitational 
energy of the galaxy, not only to heat the gas, but also to expel 
large quantities of it out of the galaxy (e.g. Bower et al. 2008). Mo- 
mentum in the relativistic jets alone is not sufficient to expel (accel- 



erate to escape velocity, ~ few o") a considerable amount of mass 
(> Mq). The maximum momentum that can be released in a rela- 
tivistic jet (which is the expected case for an outflow launched by 
BH accretion), t^Mbhc, is generally smaller than MqG, the approx- 
imate momentum needed to expel large quantities of gas ' . 

Silk & Nusser (2010) argue that radiation momentum is also 
incapable of accounting for the correlation. Let us consider King 
(2003) as an example of a model based on expulsion of gas by radi- 
ation momentum. In that model, an expression for the BH mass is 
derived (equation 15 there) where a mass of Mbh = L5 x 10 8 M:, 
is obtained for o~ = 200 km s . For this value of o", the gas mass 
inside a typical radius of 10 kpc is found to be 3 x 10 1() Mq ac- 
cording to King's equation (4). If this is the amount of gas expelled 
by radiation momentum M expe i, then the ratio of the expelled mass 
to the BH mass is M expe i/MBH ~ 200. while from This ratio from 
observations it is in fact ~ 500. This ratio is probably even much 
lower in King's model as he takes the cosmological value of 0.16 
for the baryon fraction, whereas this fraction is much larger (even 
close to unity) in the bulges of spiral galaxies. Another problem 
(though less severe) with this model is that one does not expect the 
ransfer of radiation momentum to mechanical momentum to be 100 
per cent efficient. Thus, we conclude that radiation momentum can 
expel only a small amount of gas, and that the correlation cannot 
be explained in King's model. 

On the other hand, there is more than enough energy released 
from the accretion process to expel the gas (by a factor of ~ 50), 
or rjMftnc S> MqO . The problem with energy is the opposite: if 
deposition of jets' energy into the ISM is too efficient, almost no 
star formation will occur. 

Thus, a successful theory should be based on a feedback 
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1 For a typical mass ratio of Mq ~ 10 3 Mbh. and a velocity dispersion of 
o~ ~ 200 km s ; r/ ss 0.1 is the energy fraction that is liberated by the 
accreted mass. 
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mechanism in which the deposition of the jets' energy to the ISM is 
regulated; it should also be able to derive the mathematical form of 
the SMBH-galaxy relation from the basic properties of the galaxy 
and its active galactic nucleolus (AGN). The penetrating-jet feed- 
back mechanism (Soker 2009) has these two attributes. In this 
model, the energy transfer from the jets to the ISM is efficient 
enough to expel large amounts of gas and suppress star forma- 
tion only when the jets are stopped within the bulge, and do not 
propagate to large distances by penetrating through the ISM. Soker 
(2009) found that this requirement leads to a proportionality rela- 
tion between the BH mass and the momentum parameter 

IL=M G - (1) 

c 

as detailed in Section 2. 

Motivated by the recent data and the penetrating-jet feedback 
mechanism, we carry out the present study. In Section 3 we exam- 
ine the correlation between Mbh and fi. We ignore the differences 
between elliptical galaxies, classical bulges and pseudo-bulges (e.g. 
Gadotti & Kauffmann 2009; Nowak et al. 2010; Hu 2009). In Sec- 
tion 4 we apply the results to the jet penetrating model and derive 
the approximate amount of momentum carried by the jets. We sum- 
marize in Section 5. 



2 THE PENETRATING-JET FEEDBACK 
MECHANISM 

In the penetrating-jet feedback mechanism the energy transfer from 
the jets is efficient enough to expel large amounts of gas and sup- 
press star formation only when the jets are stopped within the bulge, 
and do not propagate to large distances by penetrating into through 
the gas. The condition for jet-stopping is that the time required for 
the jet to propagate through the surrounding gas and break out of it 
must be longer than the jet crossing time (the time it takes a blob 
of material to transverse the jet's width). This situation is analo- 
gous to an attempt to pierce a hole through a plank of wood with 
a drill, when the drill moves horizontally and does not stay over 
one spot for enough time. This is referred to as the non-penetration 
condition. 

Mathematically, the jet crossing time at a typical radius r s 
(where the cooling surrounding mass M s resides) is 

2r s tan0 2r s 6 

T c = — RJ — 

V rc l 



1'iel 



(2) 



where is the half opening angle of the jet and v re ] is the rela- 
tive transverse motion between the inflowing mass and the SMBH 
( v rel ~ O). We assumes the jet to be narrow, thus tan 6 ~ 9. 

The penetration time is T p = r s /vh where vj, is the velocity of 
the head of the jet, which is derived from momentum balance (ram 
pressure balance of the jet and the ambient gas), and is given by 



I'h : 



2MfVffT 
M s 6 2 ' 



(3) 



where the subscript 'f (for 'flow') is used to denote magnitudes 
associated with the jets (thus Aff is the mass outflow rate of the two 
jets together and Vf is its velocity). M s is the mass inflow rate from 
~ 1 kpc scale. Equation (3) was derived under the assumptions of: 
supersonic motion, Vf S> vj,, and the stronger assumption that the 
inflowing speed of M s is <7. 

The non-penetration condition T c < T p leads to the following 



inequality 

. . OC 
M s> -5-87J p M B h 

v rel 



(4) 



where r\ p relates the rate of mass outflow Mf to the rate of accretion 
or BH growth Mbh 



pi = MfVf = Tj p M BH f. 



(5) 



pf is the rate of momentum carried by the two jets (that we term 
momentum discharge). Note that equation (5) is equivalent to the 
familiar relation for the total energy transfer 

2 



£bh = r]M BH c z 



(6) 



where £bh is the total energy released by the BH per unit time, and 
tj ~ 0. 1 is the accretion efficiency. The total energy is released in 
both radiation and mechanical energy of the jets. The fraction of 
kinetic energy carried by the jets is considered in Section 4. 

Equation (4) must in fact be an approximate equality: if the in- 
flow rate is above the value of the right hand side, the deposition of 
energy by the jets is efficient enough to expel the mass back to large 
distances and heat it. The inflowing mass that is not expelled by the 
jets is assumed to form stars in the bulge. Thus, time integration of 
both sides of the equation (4) (with equality sign; we also substitute 
V re ] = <T), leads to the following relation (Soker 2009, 2010) 

M «H = ^— fl, (7) 

where jl is the momentum parameter (defined in equation 1) which 
has units of mass. 

Feedback mechanisms have been discussed by many au- 
thors in the past, both to suppress gas cooling in cooling 
flow clusters (e.g. Binney & Tabor 1995; Nulsen & Fabian 2000; 
Reynolds et al. 2002; Omma & Binney 2004; Soker & Pizzolato 
2005) and to suppress star formation during galaxy formation 
(e.g. Silk&Rees 1998; Fabian 1999; King 2003; Croton et al. 
2006; Bower et al. 2008; Shabala & Alexander 2009; Soker 2009, 
2010). The penetrating-jet feedback mechanism does not make use 
of the Eddington luminosity limit, while some authors do (e.g. 
Silk&Rees 1998; King 2003). Most models (e.g. Silk&Rees 
1998; Fabian 1999) do not consider the geometry explicitly, while 
in the penetrating-jet feedback mechanism the geometry of the nar- 
row jets and the motion of their source are key issues; these intro- 
duce the factor 1/8 in equation (7). 

In the penetrating-jet feedback mechanism, the fast jets have 
two modes of interaction with the surrounding gas: if the jets pene- 
trate through the ISM gas, they deposit most of their energy at large 
distances, and thus cannot expel gas; if they cannot penetrate, they 
are shocked, and a hot bubble is formed. If the radiative cooling 
time of the hot bubble is longer than the flow time (or acceleration 
time), it efficiently accelerates the surrounding gas and expels it. 
Namely, in order to suppress star formation it is necessary that the 
jets do not penetrate through the ISM. For that, the jets should en- 
counter new material before they escape. This requires that there is 
a transverse velocity component between the jets and the ambient 
gas with which it interacts. The physical condition for efficient ex- 
pulsion of the gas is that the typical time for the jets-ISM relative 
motion to cross the jet's width in the transverse motion at a radius 
r s (from the central source) should be shorter than the penetration 
time at radius r s . 
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3 THE CORRELATION 



3.1 Methods and Sample 

The Mbh-M relation is initially assumed to be a power law of the 
form 



log 10 M BH = a±j3log I0 /i. 



(8) 



We use a least squares estimator of linear relations like f itexy of 
Press et al. (1992), that takes into account measurement errors in 
both coordinates. Like Tremaine et al. (2002), we add a constant 
residual error £o {intrinsic scatter) until Xr = 1 is obtained (see 
appendix). However, while Tremaine et al. (2002) have assumed 
the intrinsic scatter to be on the BH mass alone, we consider other 
possibilities as detailed below. Novak et al. (2006) found based on 
Monte Carlo simulations that this method estimates the slope with 
the least bias and variance. The maximum likelihood estimator used 
by Giiltekin et al. (2009) gives similar results to the above, but there 
is no freedom to manually set £y and study how it changes the re- 
sults. The ability to vary £q (or alternatively, %£) is important be- 
cause the result of the fit is sensitive to the stated measurement 
errors. First we use this method to show that the slope /3 of the re- 
lation in study is reasonably close to 1, and later we force j8 = 1 
and find the intercept a. 

We use three models for the intrinsic scatter: in the y-scatter 
model we assume (like Tremaine et al. 2002) that £o is the residual 
variance in log 10 MBH (which is traditionally the y coordinate); in 
the x-scatter model, £o is the scatter in the galactic property in ques- 
tion, in our case log 10 /x; in the orthogonal scatter model, residual 
errors are added to both coordinates, such that the combined error 
is £o in the direction perpendicular to the ridge line of the relation. 
The estimators of the three models are given in the appendix. The 
actual value of £o is meaningful only within the three models, and 
it is generally wrong to compare the values obtain from each one. 
We also note that the assumption that the scatter (in whatever direc- 
tion) is constant throughout the relation is made out of ignorance 
and may not represent the real situation. 

It is impossible to avoid specifying the direction of the intrin- 
sic scatter (Novak et al. 2006), and an extremely biased result is 
obtained if the calculation is performed under the assumptions that 
the residual variance is in the wrong coordinate. It is also claimed 
by Novak et al. (2006) that the direction of the scatter depends on 
whether it was the BH which affected the host galaxy property or 
vice versa. In the penetrating-jet feedback mechanism, the former 
is true (because the jets launched by the BH affect the gas in the 
galaxy), and thus the x-scatter model better suits the theory. How- 
ever, acknowledging that both Mbh and the momentum parameter 
are affected in a non trivial manner by phenomena such as mergers, 
we do not give the x-scatter model any preference, and regard the 
orthogonal scatter model as a compromise between the other two. 

For simplicity we use the S sample of Giiltekin et al. (2009), 
which includes 49 measured BH masses but no upper limits. In 
that paper, two different mass measurements are stated for both 
NGC1399 and NGC5128; we take the (geometric) average mass of 
each one, and the uncertainty ranges are combined. The Mq values 
are taken from table 3 in Feoli et al. (2010). 

Since our estimators cannot deal with asymmetric measure- 
ment errors, we adopt the common practice of taking 

8y = \ [log 10 M high - log 10 M low ] (9) 

where 8y is the uncertainty in the logarithmic BH mass, and M^igh 
and M\ 0VJ are the published limits of the 68 per cent confidence 
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Figure 1. The best fitting slope j3 in equation (8) as a function of the de- 
manded Xr m the three scatter models. The solid line is the x-scatter model, 
the dotted line is the orthogonal scatter model, and the dashed line is the 
y-scatter model. The two vertical lines mark the 68 per cent confidence in- 
terval of the (reduced) x~ -distribution with 47 degrees of freedom, centred 
at Xr = 1 ■ Note that this interval does not actually give the uncertainty on 
j3, which is calculated according to the algorithm of Press et al. (1992). The 
best fitting slope and intercept for Xr — 1 are given in Table 1 . 



Table 1. The best fitting values for equation (8) in the three models of in- 
trinsic scatter. In all models, £o was increased until Xr = 1 was achieved. 
The values for the intercept a vary greatly and have large errors, but all 
slopes are consistent with 1 to within two standard deviations. Note that £o 
has a somewhat different meaning in each model (see text), and does not 
indicate goodness-of-fit. 



Model 



a±8a 



13±SP 



£() 



x-scatter 
y-scatter 



0.2 ±0.6 1.04 ±0.08 0.32 
1.4±0.5 0.88±0.07 0.31 



orth-scatter 0.9 ±0.6 0.95 ±0.08 0.32 



region. The error on all Mq values is taken to be 0.18 dex; this 
dominates over the error in a (which is typically ~ 0.02 dex) in the 
uncertainty in momentum parameter. 

3.2 Results 

Fig. 1 shows the best fitting slope obtained for equation (8) using 
the three different scatter models, as a function of the achieved %*, 
The solid line is the x-scatter model (scatter in the momentum pa- 
rameter), the dotted line is the orthogonal scatter model, and the 
dashed line is the y-scatter model (scatter in the SMBH mass). In 
the point on the right where the three lines meet, x 1 is maximal and 
there is no residual variance in any direction (£q = 0). The two ver- 
tical lines mark the 68 per cent confidence interval of the (reduced) 
X 1 -distribution with 47 degrees of freedom, centred at Xr = 1- The 
best fitting values for a and |3 for %* = 1, and the obtained £o, are 
given in Table 1 . While the values of a vary greatly and have large 
errors, the errors on all slopes are at the ~ 8 per cent level, and all 
are consistent with 1 to within two standard deviations. 

This result justifies trying to modify the assumption that the 
Mbh^M relation is a power law, and instead assume that it is a linear 
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Figure 2. The Mbh~|U correlation, jl is the momentum parameter de- 
fined in equation (1). Both axes are logarithmic in M@. The data from the 
Giiltekin et al. (2009) S sample. The solid line and the dashed line are the 
best fitting straight lines according to the x- and y-scatter models respec- 
tively. The dotted line here has a slope of j3 = 1 and the intercept is given 
by equation (10). 



relation (/3 = 1). We used a similar method to f itexy to obtain a, 
as detailed in the Appendix. The forcing of the slope to 1 renders 
the three scatter models equivalent. The residual error is taken to 
be the average value obtained in Table 1 . The resulting intercept is 



a = 0.51 ±0.07. 



(10) 



The uncertainty in a is obtained using this method is much smaller 
than before, because there is much less freedom to change a while 
minimizing % ■ The SMBH mass is plotted as a function of the 
momentum parameter in Fig. 2, where the best fitting lines are 
also shown. The solid line corresponds to the x-scatter model, the 
dashed line corresponds to the y-scatter model, and the dotted line 
has a fixed slope set to j3 = 1 and an intercept given by equation 
(10). 

The motivation to check the M^n-^i correlation comes from 
the penetrating-jet feedback mechanism that gives both the slope |3 
and the intercept a from basic parameters of the systems; one of 
these parameters is determined here (see Section 4). Other models 
exist, but in many the slope and/or the intercept are not given from 
basic parameters of the system. Hopkins et al. (2007), for exam- 
ple, take the accretion rate based on an Eddington-limited prescrip- 
tion based on Bondi-Hoyle-Lyttleton accretion theory. The Bondi- 
Hoyle-Lyttleton accretion theory has the disadvantage that it can- 
not maintain a feedback in real time (Soker 2010). Their theoreti- 
cal explanation is based on pressure-driven wind, and they expect 

1 /2 9 

a relation of the form Mbh <x M G o"-. This is not compatible with 
observations, as they themselves find. 

3.3 Comparison 

We wish to compare the Mbh-M correlation we showed with corre- 
lations with other quantities such as <y. However, by setting % r — 1 
to account for the intrinsic scatter, one makes it impossible to use 
common goodness of fit tests. £o itself would be a fair indicator of 
goodness of fit only if the scatter is in the same quantity, i.e. Mbh, 
in each relation; but we have no good indication that that is the 
case. Let us nevertheless investigate other relations using the same 



Table 2. The slope (with error) and intrinsic scatter in correlations of the 
BH mass with four parameters. All relations are assumed to be linear in the 
log-log plane. The results were obtained using the y-scatter model, which 
we do not prefer. The estimated error in £o in all cases is 0.05 (see text), 
indicating that all correlations are equally good for the current data. 



Parameter 



p±Sj5 



f(i 



l" 


0.88 ±0.07 


0.31 


a 


4. 18 ±0.37 


0.38 


M G 


1.07 ±0.09 


0.33 


M G a 2 


0.74 ±0.05 


0.31 



fitting method described in Section 3.1 and the y-scatter model (so 
that the scatter is in the BH mass). We use the same dataset and 
assume all relations we test have the same form as equation (8). 

Table 2 gives the slope (with formal error) and the intrinsic 
scatter for different popular relations. The slopes were calculated 
with the y-scatter model, which we do not prefer, so that the scatter 
is in log] () A/BH in all tested correlations. We note that our results 
for the MftH-o relation are very close, but not identical, to those ap- 
pearing in table 6 of Giiltekin et al. (2009), where j5 = 4.06 ± 0.37 
and £o = 0.40 (third row there, titled T02ind). These numbers were 
obtained using the same method and data. However, it is possible 
that the discrepancy comes from different treatment of NGC1399 
and NGC5 128 (the galaxies with two Mbh measurements for each); 
Giiltekin et al. (2009) take all measurement and give each one half 
the weight, while in this work we average the values and expand 
the uncertainty ranges. We note that the sample of Giiltekin et al. 
(2009) is basically the same as that of Graham (2008b), who find 
the slope to be 4.28, which is somewhat closer to the value we find. 

All values of £o for the tested correlations are comparable; we 
give a crude estimate of the error in £o to be ~ 0.05 in all four cases. 
This is determined by using the formula 



5£o : 



3[eoteow)- £ o(Zh 2 igh)]' 



(11) 



where £o{X\ ov/ ) is the intrinsic scatter that produces a % at the low 
limit of the 68 per cent confidence range of the % -distribution 
(shown in Fig. 1); £o(Xhi 2 h) corres P°nds to the high limit of this 
range. The formal error in the scatter in the M^u-a relation as 
found by Giiltekin et al. (2009) using their maximum likelihood 
method is ~ 0.06. 

Given that the shape of the scatter is unknown and that the 
assumptions made were very simplistic, the uncertainty in £o may 
be much higher than our estimate. Thus, all tested correlations seem 
to be equally good. 



4 IMPLICATIONS FOR THE PENETRATING-JET 
FEEDBACK MECHANISM 

We now use the results of Section 3.2 to estimate the value of 
rip. We recall that r\ p is defined in equations (7) and (5), and is 
the ratio of the momentum discharge carried to Mbh c - Assum- 
ing a proportionality relation between Mbh and fl, the ratio is 
M B h/M = 10" = 3.26 ±0.51 according to equation (10). Substi- 
tuting this value in equation (7) gives 



T] p = 0.038 ±0.006. 



(12) 
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We thus derived the average value of a fundamental parameter char- 
acterizing jets launched by an AGN during the phase of SMBH 
growth at galaxy formation. 

From equations (5) and (6) it is possible to derive the fraction 
/ of the total energy that is carried by the jets. Namely, the kinetic 
power of the jets divided by the total power of the accreting SMBH 
as given by equation (6). Substituting the relativistic expressions 
for the jets' kinetic energy and momentum, we find 



(13) 



where Vf is the velocity of the jets. In the relativistic limit (Vf — > 
c), the fraction / of kinetic energy approaches a constant rj p /ri = 
0.38 (for T] = 0.1); in the non-relativistic limit (v < 0.5c) equation 
(13) becomes / = -S&— ■ Namely, in the relativistic limit and during 
the phase of SMBH growth and galaxy formation, the jets carry an 
average ~ 0.4 of the SMBH power, while in the non relativistic 
limit it is a small fraction. 

There is a lower limit set by the requirement that the jets have 
enough energy to expel the required amount of gas from the newly 
formed galaxy; the requirement is that E^f > ^MqG 2 , or 




/> 



M G a z 



(14) 



2r}M BH c 2 ' 

Using the same typical numbers as in Section 1, we find / > 0.002. 
This is an important result: even if the jets are not relativistic and 
carry only a small fraction of the energy, then this mechanism can 
work. Non-relativistic wide outflows, termed slow massive wide 
(SMW) outflows, can be formed by disk winds instead of highly 
relativistic jets from very close to the SMBH, or from a narrow rel- 
ativistic jet that very close to the SMBH turns into a non-relativistic 
jet by interacting with the ambient material (Soker 2008). 

Let us elaborate on SMW outflows. Feedback based on slow 
(~ 10 km s _1 ) massive wide jets has been applied in cooling flow 
clusters (Sternberg et al. 2007; Soker et al. 2010) and in elliptical 
galaxies (Ostriker et al. 2010) to heat and expel gas. The proper- 
ties of such SMW outflows have been recently deduced (Moe et al. 
2009; Dunn et al. 2010). For these outflows to be energetic, they 
must result from accretion close to the SMBH, rather than from 
winds, which originate in larger distances (where the potential well 
is shallow). Soker (2008) argued that fast jets launched from the 
vicinity of the SMBH can form energetic enough SMW outflows, 
and that the condition for the fast jets to power the SMW outflow is 
that the fast jets do not penetrate through the gas (as in the penetrat- 
ing jet mechanism). Namely, the mechanism studied in this paper 
to establish the Mbh~M correlation can also explain the energy of 
SMW outflow. Note that we do not confront the question of energy 
transfer from the accretion process to the SMW outflow. 

It is noted again that in the full relation between Mbh and Mq, 
the ratio of the transverse velocity to the dispersion would appear 
squared (cf. equation 4). It is assumed in the model that v re i = a. 
However, they need not be exactly equal, but we do expect them to 
be comparable. This term will surely introduce a large scatter in the 
relation; it might also introduce a systematic shift. In such a case 
the average value we derive here for r\ p incorporates this systematic 
shift. 



5 SUMMARY 

We examined the correlation between Mbh, an d the- momentum 
parameter, fl = MqG /c. The motivation for this study is the 
penetrating-jet feedback mechanism that predicts such a correlation 
(Section 2; Soker 2009), although the SMBH determines the cor- 
relation by depositing energy into the ISM rather than momentum. 
Using the sample of Giiltekin et al. (2009) we examined the corre- 
lation with jl in Section 3. We applied our results to the penetrating- 
jet feedback mechanism in Section 4. The main results of our study 
and the implied insights can be summarized as follows: 

(i) Despite the large statistical uncertainties, recent data suggest 
that the masses of SMBHs are indeed correlated with the momen- 
tum parameter of the bulge. Moreover, we find that the relation is a 
linear one. 

(ii) The penetrating-jet feedback mechanism (Soker 2009) is 
compatible with such a relation. 

(iii) In the jet-penetrating feedback mechanism, r\ p is a funda- 
mental parameter that must be obtained from observations. It has 
the same role as tj in equation (6), the relation for the total energy 
released by the accreting BH, but for the momenta of the two jets. 
We find r\ p ~ 0.04, which is of course an average value over many 
systems, and average over time in each system (during the phase of 
SMBH growth). 

(iv) This allows us to scale equation (7) by 

^l=48o( 'yUJk. 

M m V200kms- 1 / V0.04/ 



(15) 



(v) We further analysed the implication of r\ p ~ 0.04, and de- 
rived an expression for the fraction of the SMBH power that is car- 
ried by the jets (equation 13). While relativistic jets carry ~ 38 per 
cent of the released power, non relativistic jets carry a small frac- 
tion of the power (the rest is in radiation). Although the latter only 
carry a small fraction of the total power, they are still able to main- 
tain the feedback process. In this model, there is no need for "fine 
tuning" in using only a small fraction of the jets' energy, but instead 
a substantial fraction of the energy is used. 

(vi) We cannot determine with current data which of the corre- 
lations of Mbh (with /x, rr, Mq, or MqCS 2 ) is better. 
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APPENDIX 

In the fitting method we used in this paper, the following expression 
was minimized to find a and j8 : 

x ,e ie 2+j3 2 4+E 2 (/Mo)' 

where Xj and y,- are the data points, £„■ and £„■ are the uncertainty 
values for galaxy i; E (|3, £o) is different in the three scatter model 
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used 



E 2 (j3,eb) ; 



2 C 2 



Fti 



](1+P 4 )/(1+P 2 



x-scatter 
y-scatter 
orthogonal scatter. 



(2) 



If no intrinsic scatter is considered, then E = and the % esti- 
mator is symmetric in x and y, as noted by Tremaine et al. (2002). 
When adding the intrinsic scatter, only the orthogonal model pre- 
serves this feature. In the results shown in Fig. 1 , we increase (from 
zero) the value of £q so that % P er degree of freedom (denoted by 
X?) approaches its expectation value of unity. When demanding 
j8 = 1 (as we do in obtaining equation 10) the three models become 
identical. The calculation of the formal error on the best fitting pa- 
rameters is based on finding an (a, j3) pair for which x 2 is larger by 
1 from its minimal value, as detailed in great length in Press et al. 
(1992). 

One can think of orthogonally scattered data in the following 
way: a certain law of nature produces a linear relation between x 
and y, but various other natural processes cause objects to slightly 
move on the xy plane. Points (x,-,y,) with no measurement errors 
will have distances d\ from the ridge line of the relation, which 
are distributed normally with variance £q . Note that while the units 
of £o are understandable in the other models (as they correspond 
to one axis), here the units of the intrinsic scatter are intermediate 
between the two coordinates. Thus, the obtained value of £o is only 
meaningful in a particular plane, and cannot be used to compare 
goodness of fit of different relations. 

We did not study the orthogonal estimator thoroughly, and it is 
proposed here just as a logical compromise between the two other 
choices, we briefly describe the behaviour of this estimator: when 
the slope is high, the orthogonal scatter model behaves like the x- 
scatter model, and in fact one usually gets a slightly higher best 
fitting /3; when the slope is too low, then the estimated /3 is close 
to that obtain by the y-scatter model. Therefore, it seems that this 
approach is good primarily when the slope is moderate, which was 
the case studied here. 
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